ABSTRACT
In this paper, we will use the ordinary differential equations for the SIR model as a non-linear system with the Runge-Kutta numerical method of the 6th order and 7th order to generate default values for (susceptible people, infected people and recovered from disease) for epidemic disease COVID-19 by giving it the initial values for a population in a particular country. The MATLAB program was used to solve the two problems (6th & 7th order) and obtain the results.Through this work we notice the difference between the results of the two methods and the solution period as well as the estimated error value for the solution in each problem and the comparison between results and solutions for both methods and shown in a table for clarity.We also used the binary test (0-1) to know the behavior of the disease in terms of chaos and the results indicated that this disease(covid-19) is chaotic and irregular and by using the MATLAB program we obtained the figures and results and that illustrate its chaotic behavior (kcorr 6th =0.9212) and (kcorr 7th = 0.9560),as well as, the figures that illustrate the system's workflow in this paper. © 2023 American Institute of Physics Inc.. All rights reserved.